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Abstract 

In the framework of the Constrained Minimal Supersymmetric Standard Model (CMSSM) 
we evaluate the electroweak fine tuning measure that provides a quantitative test of super- 
symmetry as a solution to the hierarchy problem. Taking account of current experimental 
constraints we compute the fine tuning at two-loop order and determine the limits on 
the CMSSM parameter space and the measurements at the LHC most relevant in cover- 
ing it. Without imposing the LEPII bound on the Higgs mass, it is shown that the fine 
tuning computed at two-loop has a minimum A = 8.8 corresponding to a Higgs mass 
rrih = 114 ±2 GeV. Adding the constraint that the SUSY dark matter relic density should 
be within present bounds we find A = 15 corresponding to rrih = 114.7 ± 2 GeV and this 
rises to A = 17.8 (rrih = 115.9 ±2 GeV) for SUSY dark matter abundance within 3er of the 
WMAP constraint. We extend the analysis to include the contribution of dark matter fine 
tuning. In this case the overall fine tuning and Higgs mass are only marginally larger for 
the case SUSY dark matter is subdominant and rises to A = 28.7 (rrih — 116.98 ± 2 GeV) 
for the case of SUSY dark matter saturates the WMAP bound. For a Higgs mass above 
these values, fine tuning rises exponentially fast. The CMSSM spectrum that corresponds 
to minimal fine tuning is computed and provides a benchmark for future searches. It is 
characterised by heavy squarks and sleptons and light neutralinos, charginos and gluinos. 
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1 Introduction 

With the Large Hadron Collider up and running, the search for TeV-scale SUSY is now 
significantly closer. Many models of physics beyond the Standard Model (SM), in particular 
its minimal supersymmetric version (MSSM), will be directly tested. To do so, one has to 
quantify the viable range of parameters entering in these models, which impact in particular on 
the scale of its low-energy supersymmetric spectrum. Previous theoretical and experimental 
constraints such as those from LEP have already tested a considerable amount of the MSSM 
parameter space and identified bounds on it. To study these bounds further it is instructive 
to investigate from a quantitative perspective the impact of the, so far, negative searches for 
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low-energy supersymmetry. A quantitative measure of this impact is the fine-tuning measure 
A [TJ [2] , that quantifies the degree of cancellation between unrelated parameters that is needed 
to fix the electroweak scale and can be extended to include the fine tuning needed to obtain 
an acceptable dark matter abundance. In this paper we shall perform such an investigation, 
computing A at two- loop leading log order in the Constrained MSSM (CMSSM). For previous 
studies of the fine-tuning problem in a similar context see [T]-|18j. 

The electroweak fine-tuning measure, A, provides a measure of the probability of unnatural 
cancellations of soft masses in the expression of the electroweak scale v 2 ~ — ^ m 1oft «M> 
is Higgs quartic coupling) after including quantum corrections. So A measures the stability 
of the MSSM electroweak scale at the quantum level, with all available experimental and 
theoretical constraints imposed. These include the LEP mass bounds on supersymmetry 
masses, charge/colour breaking constraints, the dark matter relic density constraint, and the 
measurement of b — > 57 and the anomalous magnetic moment of the muon. In what follows we 
identify the constraints with the largest impact on A. We also extend the analysis to include 
the fine tuning, A^, needed to satisfy the constraints on the SUSY dark matter abundance. 
The method can be readily extended to other models that claim to solve the hierarchy problem. 
For the case that the fine tuning is reduced by new states with mass well above the TeV range, 
one may extend the analysis using the effective Lagrangian in which the very heavy states 
have been integrated out. For fine tuning and related issues in such scenarios see for example 
[17 } I20 | \21 \ \22 \ I23j . where the MSSM Higgs mass can be increased nearer the LEP bounds by 
classical effects due to new physics beyond few TeV, which ultimately reduces the fine tuning. 

The fine-tuning problem in the MSSM is important not only for supersymmetry searches, 
but also for the Higgs physics since it is intrinsically related to the value of the lightest Higgs 
mass rrih, currently restricted by the LEPII lower bound of 114.4 GeV |19j . As a result 
searches for rrih are relevant to supersymmetry phenomenology. In particular, the need to 
increase the SUSY prediction for the Higgs mass by radiative corrections above the LEPII 
bound means that the electroweak fine tuning measure rises exponentially with the Higgs 
mass. If A becomes too large one can conclude that SUSY fails to provide a solution to 
the hierarchy problem. The parameter configuration (consistent with the non-observation of 
SUSY states) that minimises A gives an indication of its most likely values. We identify this 
configuration and investigate its phenomenological implications. Also for a given upper value 
of A one can extract the corresponding range of parameter space of the CMSSM and of the 
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superpartners masses. 

In Section 2 we present the calculation of the electroweak fine tuning measure to two loop 
order. We are not aware of a similar analysis of the fine tuning problem at this level of precision 
(two- loop leading log), which is responsible for a A smaller than that usually found in the 
literature. In Section 3 we discuss the dependence of the fine tuning on tan/? and the impact 
on the fine tuning coming from imposing the bounds on the SUSY spectrum and the limits on 
b —7- sj. Using the dependence of the fine tuning measure on the parameters of the CMSSM 
we then determine their allowed range consistent with a given value of A. This provides a 
quantitative measure of the remaining parameter space range that remains to be tested. Next 
we determine the dark matter abundance as a function of the fine tuning measure showing 
that low fine tuning is consistent with acceptable SUSY dark matter abundence. We conclude 
this Section by considering the implications for the Higgs mass following from requiring low 
fine tuning. We show that, without imposing the LEPII bound on the Higgs mass, A has a 
minimum for a region of nih near the LEPII bound. In Section 4 we extend the fine tuning 
analysis to include the fine tuning needed either to satisfy the dark matter bound or to 
saturate the bound with SUSY dark matter. Finally in Section 5 we discuss the predictions 
for the superpartner mass spectrum from the fine tuning bound. We also determine the most 
likely spectrum that minimises fine tuning and discuss the relative importance of various LHC 
measurements in the test of the CMSSM. Section 6 presents a summary and our conclusions. 

2 Computing electroweak scale fine-tuning A at two-loop level 

In this section we present the strategy for evaluating the MSSM fine-tuning at tree, one-loop 
and two-loop (leading log) level; particular attention is paid to clarifying the impact on fine 
tuning of the quantum corrections to couplings and masses. With the standard two-Higgs 
doublet notation, the scalar potential is 

V = m\ \Hx\ 2 +m| \H 2 \ 2 - (m§ H x ■ H 2 + h.c.) 

+ — Ai |i^i| 4 + — A 2 I H 2 1 4 + A3 I -Hi 1 2 I -£^2 1 2 + A4 \ H\ ■ H 2 \ 2 



+ 



X - A 5 (#! • H 2 ) 2 + A 6 |#i| 2 (tfi • H 2 ) + A 7 \H 2 \ 2 ■ H 2 ) + h.c. 



(1) 

The couplings A,- and the soft masses receive one- and two-loop corrections that for the MSSM 
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are found in [24|. [25] . We shall use these results to evaluate the overall amount of fine-tuning 
of the electroweak scale. Technical details of the procedure can be found in the Appendix. 
To evaluate the fine-tuning, it is convenient to introduce the notation 

m 2 = m\ cos 2 P + m\ sin 2 P — m\ sin 2/3 

h C os 4 /3+^ sin 4 /3 + ^- 5 ~- 2 ™ ' ™o*f* —2* ■ * ^ 



A = -1 cos 4 p + -f sin 4 p + -2S sin 2 2/3 + sin 2/3 (A 6 cos 2 p + A 7 sin 2 p) (2) 

IrJJ tl 



with the assumption that, in the MSSM, at the UV scale m\ = m\ = m\ + Uq while m 2 = 
BofiQ. The couplings A- are assumed to be real and A 345 = A 3 + A 4 + A 5 . the MSSM, at 
the tree level they are 



A 1 = A 2 = l/4( 5l 2 + <? 2 ), A3 = l/4( 5 2 - 5l 2 ), A 4 = -l/2^, A 5 , 6 ,7 = (3) 

where are the C/(l) and SU(2) gauge couplings respectively. 

The fine-tuning amount wrt to a set of parameters {p} of the MSSM is then [TJ [2] 

<91nf 2 

A = max |A p | p={ ^ mg m?/2iA g A p = (4) 



where all p are input parameters at the UV scale of CMSSM, in the standard notatioio The 
minimisation of V gives 

v = -m /A, 2A^ =m- (5) 

which fix v and P as functions of the above MSSM bare parameters. Taking into account that 
m 2 = m 2 (p,P), A = X(p, P) we can find dp /dp from the second minimum condition for V. 
Using this, one findj^ [IT], see also [17] : 



P 

z 



d 2 m 2 2 d 2 \\ f d\ 1 dm 2 \ dm 2 d 2 \ d\d 2 m 21 



2 dp 2 +V dp 2 J \dp + v 2 dp ) + dp dpdp dp dpdp 



,d 2 m 2 2 d 2 X\ v 2 (d\\ 2 



1 When using the Yukawa couplings at the low energy scale as an input, expressed in terms of the Higgs 
vev and quark masses, then m and A pick up additional, implicit dependence on v = (h) and tan /?. Neglecting 
this dependence when evaluating A p does not bring significant changes to the results for final A. 

2 One could also include in the set of parameters p, the top Yukawa coupling or the strong coupling 013. 
For such parameters which are measured, one can use the modified fine tuning definition [26] and with this, it 
turns out that in the cases we consider their associated fine tuning never dominates. 

3 Later the min condition (fixing /3) is used to replace Bo by tan ft as an independent parameter. 
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This result takes into account the dependence of /3 on the MSSM set of parameters (p). This 
formula also takes into account the loop-effects to the quartic couplings as well as the tan /3 
dependence of the radiative corrections on the parameter "p". As we shall see later, such 
effects tend to reduce fine-tuning, in many cases rather significantl}0. 

For comparison with similar studies of fine-tuning, a comment is in place here. After some 
algebra, one can show that the general formula of A (eqs.(H|), ©) reduces, in the limit of 
removing the loop corrections to quartic couplings Aj, i = 1,2..., to the more familiar "master 
formula" [28] (see also [HE]) 



p f dm\ 2 dm 2 
A„ = 5— — — — 2" i —5. tan P ■ 



(tan 2 13 — 1) m z [dp dp 



tan B \ m 2 7 
1 + 



cos 2/3 



m 2 + rh 2 



dm 2 - „ (dm 2 dmk 

3 sin 2/3 — -i + 



dp V dp dp 



(7) 



This formula is sometimes used as the starting point in works that evaluate electroweak scale 
fine-tuning, by using in it the loop corrected soft masses. However, for accurate estimates, 
it is necessary to take full account of radiative corrections, as done by eqs.(j3|), ©. Indeed, 
the loop corrections to the quartic couplings significantly reduce the amount of fine-tuning, 
in some cases by a factor as large as 2, and these corrections are not accounted for by eq. ((TJ) . 
This can be seen by considering the one-loop correction 5 to A2 — > A2(l + 6), due to stop/top 
Yukawa couplings. Usually 5 = 0(1). Including it one finds (for details see eq.(26) in [17]): 

Ap ^ (l + £)m 2 +°( 1 /tan/3), p = fi 2 , m 2 , m\ /2 , A 2 , £? 2 . (8) 

This is showing that one-loop corrections to the quartic coupling reduce the amount of fine- 
tuning significantly. In fact it is the smallness of the quartic Higgs coupling (fixed in the 
MSSM by gauge interactions) that is at the origin of substantial tree-level fine tuning. That 
this is so can be seen from the relation v 2 ~ —mf/X where v is of 0(100) GeV, rrii ~ 0(TeV) 
while at the same time A < 1 , which makes it difficult to separate the EW and SUSY breaking 



4 The radiative corrections to couplings \j and soft masses rrii bring about additional field dependence, and 
therefore additional v and j3 dependence. If we include the extra v dependence, we find the fine-tuning A p 
changes into A p — A p /(1 — A v 2) with A„2 defined by eq.©. |A„a| <C 1 in most cases examined. We do not 
include this effect here and work with A defined by eqs.Q, ©. 
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scales (for a discussion see [7]). Loop corrections increase the quartic couplings and in most 
cases reduce the overall amount of fine tuning. 

For a general two-Higgs doublet model formulae (jl]), © can be expressed in terms of 
only derivatives of couplings and of masses wrt to the corresponding parameter, see Appendix 
in [T7] (also Appendix A.l). For the CMSSM we use this result, in which we consider the 
full two-loop (leading log) corrections to the quartic couplings and masses. This defines 
unambiguously our procedure for evaluating the EW fine-tuning in CMSSM. 

3 Electroweak fine tuning and its effects on the Higgs mass 

In the following, the numerical results we present for A include two-loop corrections with: 

• radiative electroweak breaking (EWSB), 

• non-tachyonic SUSY particle masses (avoiding colour and charge breaking (CCB) vacua). 

• experimental constraints considered: bounds on superpartner masses, electroweak precision 
data, b — > sj, b — > fifi and anomalous magnetic moment <5a M , as detailed in Table [TJ 

• consistency of with the LEPII bound (114.4 GeV) and/or consistency with thermal relic 
density constraint, only if stated explicitly. 



Constraint 


Reference 


SUSY particle masses 
Sa^ < 366 x KT 11 
3.20 < 10 4 Br(6 -> sj) < 3.84 
Br(6 -> fifi) < 1.8 x 10~ 8 
-0.0007 <5p< 0.0012 


Routine in MicrOmegas 2.2, "MSSM/masslim.c" 
PDG (sys. and stat. la errors added linearly) 
PDG (sys. and stat. la errors added linearly) 

Particle Data Group 

Particle Data Group 



Table 1: Constraints tested using MicrOMEGAs 2.2 with SuSpect 2.41 spectrum calculator. Particle 
Data Group: |http://pdg. lbl.gov/ 



Using these constraints we evaluated A numerically. The LEPII bound on the mass of the 
Higgs provides an important constraint for the MSSM since it requires quantum corrections 
in order to be satisfied. Large quantum corrections need in turn large soft masses, which in 
turn trigger large fine-tuning. This is seen from the loop corrections to which give a strong 
exponential dependence on m^, A ~ ~ exp(m^/m^ op ). To examine this dependence in 
detail, in the following we choose to present the numerical results as a function of m^. Unless 
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Figure 1: Fine tuning vs Higgs mass, in a two-loop analysis. The data points are for 2 < tan/3 < 55. 
The solid line is the minimum fine tuning with central values («3,mt) = (0.1176, 173.1 GeV). The 
dashed line corresponds to (a 3 ,m t ) = (0.1156, 174.4GeV) and the dotted line to (0.1196, 171.8 GeV), 
to account for la experimental errors in a 3 and top mass [34 . This is the "worst" case scenario, when 
such deviations combine such as to give the largest variation of A. An increase of a^(rnz) or reduction 
of m t (mz) by la have similar effects, which can be also understood from the relation between the mass 
of top evaluated at mz and at nit- Keeping either a 3 or m t fixed to its central value and varying the 
other within la brings a curve situated half-way between the continuous line and the corresponding 
dashed or dotted line. The LEPII bound of 114.4 GeV is indicated by a vertical line. Note the steep 
(« exponential) increase of A on both sides of its minimum value situated near the LEPII bound. 

stated otherwise, the LEPII bound on is not imposed. The relic density constraint is 
imposed only after all the constraints other than the LEPII bound on m/j are satisfied and, 
when done, this is stated explicitly. 

Before proceeding to present our numerical results obtained with the above constraints, 
let us mention the details of the procedure followed. First the fine-tuning is evaluated at two- 
loop order, including the dominant third generation supersymmetric threshold effects to the 
scalar potential. The scan is done over all parameter space using a slightly simplified two-loop 
calculation performed by a Mathematica code, based on the formulae in the Appendix, and 
optimised to run quickly. For the points in phase space that have the smallest fine-tuning 
(say with A < 1000), the analysis is re-done using (the slower) SOFTSUSY 3.0.10 [29] that 
includes all the two-loop radiative effects mentioned. This two-step procedure is extremely 
important, since otherwise the CPU run time using SOFTSUSY alone would be about 6 years 
(when run on 30 parallel processors at 3GHz each), which prevented previous investigations 
at this precision level. Our two-loop analysis is also important because there is a significant 
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Figure 2: The plot displays the various contributions max|A p |, p — fig, uiq, Aq, Bq, m,y 2 , to the 
electroweak fine-tuning A presented in Figure [T] The largest of these for all mh gives the curve 
presented in Figure [T] At low m^, A^2 (red) is dominant, while at large m/,, A m 2 (blue) is dominant 
(with A^2 reaching similar values near 120 GeV). The transition between the two regions is happening 
at about 114.5 GeV. Note that in this plot the LEPII bound is not imposed at any time. Although A h 2 
(purple) is presented above for illustration, this contribution is always sub-dominant when assuming 
the modified definition of fine-tuning [26], appropriate for measured parameters (as we do in the text). 

difference between one-loop and two-loop values for overall A, and was not performed in the 
past. QCD effects can compete at two-loop with Yukawa couplings effects, can dominate 
them and also displace the minimum of the fine tuning from its one-loop value. Regarding 
the Higgs mass, its value computed with SOFTSUSY agrees with that found using SuSpect 
[30j within 0.1 GeV, but can differ by ±2 GeV [3T] from the value found using FeynHiggs [32J. 
We use the SOFTSUSY Higgs mass for all figures in the paper. Given the small discrepancy 
with FeynHiggs, coming from higher order terms in the perturbative expansion, the LEPII 
bound should be interpreted as > 114.4 ± 2 GeV. In the following analysis the results are 
always quoted with respect to the central value of m^. 

Turning now to the numerical results, Figure [1] presents the two- loop result for the depen- 
dence of overall electroweak fine tuning A as a function of the Higgs mass. The dark matter 
constraint and the LEPII bound on mh are not included. The loop effects reduce the fine 
tuning amount; the dominant effects come from quantum corrections to the quartic couplings 
(rather than to soft masses), which are increased by radiative effects and thus reduce A. A 
is seen to have a minimum close to the LEPII bound of m/,. The individual contributions 
to A are shown in Figure [2J Below the LEPII bound, detailed calculations show that the 
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minimal value of A is dominated by A^2 and this increases rapidly for decreasing mh- For 



values of mn above the LEPII bound, A is dominated bj^l A m 2. This happens at the edge of 
the focus point region. The transition from the dominant A^2 regime to the dominant A m 2 
regime occurs near the LEPII bound value, and this is the point where the QCD radiative 
effects become important. This can be seen from Figure [T] where an increase by la of 
corresponds to a larger A for a same, fixed value of m^. In this sense one could even say that 
the minimal value of A is situated at the transition region between dominant effects, Yukawa 
versus QCD interactions. Away from the minimum of A, fine tuning increases dramatically, 
roughly exponentially. This is because, as discussed above, A depends exponentially on mh- 
In conclusion, the fine tuning at two-loop with all the latest constraints is minimised for 

A « 8.8, m h = 114 ± 2 GeV (9) 

with the theoretical uncertainty of ±2 GeV, explained earlier. 

Note that our analysis also investigated the contribution to A coming from the uncertainty 
in measured "parameters" such as top Yukawa and strong coupling. Using the modifieqj 
definition [26], appropriate for measured parameters, we find their fine-tuning is sub-dominant. 

Finally, let us also mention that the reduction of fine-tuning that we have seen is mostly 
due to (two-) loop corrections considered, particularly to quartic couplings and is actually 
very significant, given the conservative scenario considered here assuming universal gaugino 
mass structure of the CMSSM; relaxing this condition could reduce |6j A further. 

3.1 Constraints on A from fixing tan/5 

It is interesting to examine the fine-tuning for fixed values of some of the parameters present, 
to see the individual impact of such constraints on A. Here we do this for a fixed value of 
tan p. This is showrj^ in Figure [3] for increasing vales of tan j3, with tan (3 = 2 (blue), tan /3 = 3 
(red), tan/3 = 10 (green). Increasing tan/3 shifts the curves of A towards larger mh and, to a 
limited extent, to lower fine-tuning, reached for medium tan f3 ~ 10, also for larger tan (3 ~ 40 
(see later, FigJ7](a)). As may be seen from the figure, the two-loop expressions for soft masses 



Larger rtih requires larger m t - ~ m , and a larger mo, above the focus point region, increases A m 2. 

6 A P = A p x (up/p) where a p is the la error in the parameter p derived from experimental observation. 

7 In Figures [3] [4] [5] only, A is computed with our Mathematica code instead of SOFTSUSY, due to long 
CPU time constraints. This explains the small difference in shape between the two-loop line in these three 
figures from the more accurate one in Figure [T] and all other figures (based on SOFTSUSY). 



9 



70 80 90 100 110 120 130 

Higgs Mass /GeV 

Figure 3: Minimum fine tuning versus Higgs mass at two-loop, showing the influence of loop effects 
and tan/3 on fine tuning. All constraints listed in Table [1] are included. The upper and lower lines 
associated with the coloured regions are the 1-loop without thresholds for A and soft masses and "full" 
1-loop results respectively (similarly for the grey region, for all tan (3). The minimum 2-loop fine tuning 
is found between these two cases. The solid lines refer to the scan 2 < tan/3 < 55. The black points 
give the positions of minimum A for fixed tan/3 from 2 to 4 inclusive in steps of 0.5. 

and couplings bring values of A which are situated between the higher "tree-level" curve (i.e. 
tree- level for Aj, one- loop for soft masses without field dependent threshold effects) and the 
lower one-loop case (one-loop for both Aj and soft masses). This is consistent with what 
one would expect from a convergent, perturbative, loop expansion. The roughly exponential 
behaviour of A at large for fixed tan/3 can bring a significant variation of A. 

3.2 Constraints on A from the SUSY spectrum 

Here we examine the impact on A due to constraints related to the supersymmetric spectrum. 
The key features of the impact of this spectrum on A can be seen from the limits on the 
chargino mass considered in Figure [H Currently, the (lightest) chargino mass bound is the 
most important, followed by that of the neutralino. It also turns out that the gluino mass 
limit is not very constraining. These results follow recent experimental data, since using the 
1998 data it was the neutralino mass bound that was more constraining for A. 

In Figure HI the effect of the chargino mass m xi > 94 and m xi > 80 GeV is shown by 
the upper and lower continuous curves. While these can have some impact on fine-tuning for 
values of mh already ruled out experimentally, for mh > 114.4 GeV, the effect is overlapping 
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Figure 4: Minimum fine tuning vs Higgs mass, showing the influence of mass constraints on fine tuning. 
The results are at 2-loop with the upper shaded (coloured) areas connecting the case of only applying 
the SUSY spectrum constraints (lower line) to that with all constraints listed in Table [T] (upper line) . 
The lower shaded (coloured) areas connect the cases of only applying a chargino lower mass limit of 80 
and 94 GeV for the lower and upper lines respectively. The results for the scan 2 < tan j3 < 55 are also 
shown by the grey shaded area, with similar convention for upper/lower continuous lines delimiting it. 

that from 6 — > sj (see later). Further, the close vicinity of the upper and lower dotted curves 
corresponding to all constraints in Tableland to the SUSY spectrum limits respectively shows 
that the latter are the main constraints at the moment for A at low tan/3. The graph presents 
A and rah computed using quartic couplings and mass expressions evaluated at two-loop. 

3.3 Constraints on A from b — > 

Figure O gives the impact of the 6—7-57 constraint on A. The lower limit of the 6—7-57 
constraint for a given coloured area (fixed tan/3) restricts the right hand edges of the plot, 
while the upper limit restricts its left hand side. These curves also depend on the mass limits - 
these constraints are not fully independent. For the experimentally allowed area of > 114.4 
GeV, the impact of the constraint 6 — > sj is rather small; in this case its effect is overlapping 
that of the SUSY mass limits, as can be seen from the rhs of the plots for individual tan/3 
plots. The combination of the SUSY mass limits and 6—7-57 constraint currently dominate 
the restriction on how small the fine tuning could be, see also Figure [6j In this last figure one 
can easily see, at two-loop, the impact on A of removing the 6 — > 57 constraint. For a related 
analysis of 6 — > S7 see recent [27] . 
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Figure 5: Minimum fine tuning vs Higgs mass, showing the influence of the b — > constraint. 
The results are at 2-loop with the lighter shading connecting the case of only applying the SUSY 
spectrum constraints (lower line) to that with also the b — > sj constraint listed in Table [TJ (upper 
line). This upper edge of this shading is indistinguishable from the solid line which includes all 
constraints in Table [TJ The darker shading extends up to the minimum fine tuning limits for the case, 
3.52 < 10 4 Br(6 — > sj) < 3.77, with the other constraints as given in Table [TJ The results for the scan 
2 < tan p < 55 are also shown, between the two continuous and almost parallel lower curves. 



The other constraints listed in Table [TJ do eliminate further mSUGRA points, but have a 
negligible effect on the fine tuning limits. With the current mass limits, a change in the 8a ^ 
constraint by factors of 2 or more does not affect these results significantly. 

3.4 Constraints on A and the CMSSM parameters 

The fine tuning measure can be easily applied to establish the remaining allowed range for 
the MSSM SUSY parameters. In Figure [7] we plotted the dependence of the total fine tuning 
wrt various parameters. To understand some aspects of the dependence of the electroweak 
scale fine tuning on the MSSM parameters, let us use, for the sake of discussion, one of 
the two minimum conditions which, when ignoring quantum corrections to quartic couplings, 
simplifies to: 

m 2 z m\ — m| tan 2 /3 



2 tan 2 p - 1 

Using 2-loop RGE solutions at tan/3 = 10 (for details see the Appendix), one has 



(10) 
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Figure 6: Minimum fine tuning vs Higgs mass, showing the influence of the b —> sj constraint. In 
the presence of this constraint with 3<r limits, the red (lower) region is removed. Within the la limits, 
the green (middle) band is removed leaving the blue (upper) points. This is a two-loop leading log 
approximation, obtained using SOFTSUSY. 



m\{m z ) ps 0.99 /x§ + 0.946 m\ + 0.331 m\, 2 +0.044 A m 1/2 -0.013^ 



ml{m z ) ps 0.99/^ -0.080 ml -2.865 m\, 2 + 0.445 A m 1/2 -0.099^ 



(11) 
(12) 



It is the large cancellation between the //q and m\, 2 terms that is often responsible for the 
large fine tuning (note however that this argument ignores the impact of quantum corrections 
to quartic couplings, known to reduce the fine-tuning). This leads to the approximate relation 



A^ 2 ~ ^m 2 ' ^ s ^ ow ^ ne tuning prefers small fi , small m l i 2 is also preferred and this is 



observed in Fig [7] (e), (f). The rise in fine tuning at small my^ is a consequence of the 
constraints such as the chargino mass limit. 

The near flat distribution of minimum fine tuning in m is a result of the coefficient of m 
in m 2 being driven close to zero. The fine tuning with respect to m then rarely dominates, 
until we reach values of above the LEPII bound (tjiq at the edge of focus point region). 
The result of applying the Higgs mass constraint also excludes a region with small ra\j 2 at 
m-Q below 1.5 TeV. The focus point at m ~ 1.5 TeV where the minimum of ra\j 2 is possible, 
corresponds to the point where fine tuning is minimised. This only occurs for large tan f3, and 
this available "dip" in fine tuning in the mSUGRA space disappears as tan/3 is reduced. 

Figure (7]^c) indicates that a small trilinear coupling | ^4 1 < 1 TeV is preferred for the 
smallest fine tuning. This follows from a similar argument for preferring small mi/2- Increasing 
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Figure 7: Dependence of minimum fine tuning on SUSY parameters (fi > 0, relic density unrestricted). 
The solid, dashed and dotted lines are as explained in Fig[TJ No bound on m,h is applied in figure (a). 
In (c), (d), (e), the darker shaded regions are eliminated when m h > 1 14.4 GeV is applied for the case 
with the central (a 3 ,m t ) values. In (b) and (f), m/j > 114.4 GeV is applied, and the points in (f) are 
only for the central (a 3 ,m t ) values. 



|^4q| requires larger cancellations with fi to set the electroweak scale. However, once the Higgs 
mass constraint is applied, A Q is driven negative for small tan f3 in order to maximise the 
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stop mixing. The related increase in the minimum fine tuning from being in this region of 
parameter space then follows. This is important for small tan j3 where the tree level Higgs mass 
is smallest. The sign structure of the UV parameter coefficients in m 2 leads to a preference 
in a small, positive A Q . 

As mentioned earlier, the fine tuning measure A can be used to establish the remaining 
parameter space of the CMSSM compatible with a solution to the hierarchy problem. As- 
suming that A = 100 is the upper limit beyond which we consider that SUSY failed to solve 
the hierarchy problem, we obtain the following bounds: 

m h < 121 GeV 5.5 < tan/3 < 55 

H < 680 GeV 120 GeV < m 1/2 < 720 GeV (13) 

to < 3.2 TeV -2.0 TeV < A < 2.5 TeV 

These values can be easily re-calculated for a different value of A. 

3.5 Constraints on A from the relic density 

It is interesting to see the impact on A and on the CMSSM parameters from the presence 
of the dark matter relic density constraint (examined using micrOMEGAs2.2 |35j ) and the 
LEPII constraint on m^. These are rather strong constraints, particularly in the "restrictive" 
context of CMSSM that we study, with universal gaugino mass. The results are presented in 
Figure where the relative impact of the LEPII constraint can be seen by comparing the 
left and right plots. If the observed dark matter abundance is imposed as a constraint on the 
CMSSM, then the range of values given in (|13p and valid for A < 100 is further restricted, as 
seen in Figure The condition that the SUSY LSP should provide the observed dark matter 
abundance as well as the constraint mh > 114.4 GeV removes the intermediate values of to^ 
and mo, but has a rather small impact on Aq. 

3.6 Prediction for m/, from minimising A and saturating the relic density 

The relic density constraint can be combined with that of minimal electroweak fine-tuning 
A to make an interesting prediction for m^. Figure shows the impact of non-baryonic 
relic density constraint on A presented in Figure [TJ Obviously, not all initial points in A 
satisfy this constraint, and this is shown in Figure [9] by the red and blue points which do not 
fill the whole area above the continuous line of minimal A. As expected, the additional dark 
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Figure 8: Dependence of minimum fine tuning on SUSY parameters, with the relic density saturated 
within 3ct of the WMAP bound (in red). The WMAP bound is tlti 2 = 0.1099 ± 0.0062 [33]. The 
blue (darker) points do not saturate the relic density fth 2 < 0.0913 (3a deviation). The impact of 
the constraint > 114.4 GeV is also considered. Compare this figure to Figure [7] where relic density 
constraint was not included. The parameters values quoted in eq. (fT3j) are further restricted, as seen 
from these plots. The continuous line is that of minimal electroweak A (no relic density constraint). 
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Figure 9: Fine timing vs Higgs mass with the influence of the WMAP bound. The minimum fine 
tuning at 2-loop for 2 < tan ft < 55 is given by the solid line when including all the constraints listed in 
Table [TJ Left figure: The blue (darker) points sub-saturate the relic density. The red (lighter) points 
give a relic density within the la bounds, tth 2 = 0.1099 ± 0.0062. The 'strips' of points at low Higgs 
mass appear due to taking steps of 0.5 in tan ft below 10. A denser scan is expected to fill in this region. 
Similarly, more relic density saturating points are expected to cover the wedge of sub-saturating points 
at m h — 114 GeV and A > 30. Right: as for left, within 3tr WMAP bound (in red). The continuous 
line is that of minimal electroweak A without the relic density constraint. 

matter constraint prefers in some cases larger A relative to its minimal value (continuous line) 
obtained only with the constraints in Table [TJ However, as can be seen in the plots, the region 
of m/j values where this constraint is indeed relevant is actually ruled out by LEPII bound 
rrih > 114.4 GeV; above this value the two curves on the boundary are almost overlapped 
and the constraints in Table [J are sufficient to also satisfy the thermal relic density; note 
that the red points in the left (right) plots in Figure [9] satisfy the relic density within la 
and 3a WMAP bounds [33], respectively. The results are obtained as usual, using two-loop 
values for the quartic couplings and soft masses and corresponding threshold effects, (with 
the SOFTSUSY and micrOMEGAs codes). 

It is important to notice that, without imposing the LEPII bound, at the two-loop level, 
the smallest fine tuning A consistent with the relic density WMAP bounds |33] predicts a 
mass for the lightest Higgs as follows: 



m h = 


114.7 ±2 


GeV, 


A = 


15.0, 


(sub-saturating the WMAP bound) 


m h = 


116.0 ±2 


GeV, 


A = 


19.1, 


(saturating the WMAP within la). 


m h = 


115.9 ±2 


GeV, 


A = 


17.8, 


(saturating the WMAP within 3a). 
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Figure 10: Left (Right) figure: Relic density fine tuning, A in function of the Higgs mass, at two-loop 
level, for a ler (3a) WMAP bound (red), respectively. The blue (darker) points sub-saturate the dark 
matter relic density. 

To conclude, minimising the fine-tuning together with the constraints from precision elec- 
troweak data, the bounds on SUSY masses and the requirement of the observed dark matter 
abundance lead to a prediction for m/j, without imposing the LEPII bound. This is an in- 
teresting result, and represents our prediction for the CMSSM lightest Higgs mass based on 
assuming A as a quantitative test of SUSY as a solution to the hierarchy problem. 



4 Dark matter fine tuning and its effect on the Higgs mass 

The dark matter abundance can be very sensitive to the choice of parameters and can introduce 
a new fine tuning to the model. To quantify this it is interesting to consider the dark matter 
fine tuning A^ wrt the CMSSM parameters, and to determine its impact on the overall fine 
tuning (for earlier studies see [15^ [TB] and references therein) . Its definition is similar to that 
of A: 

(15) 

q=mo ,tan 



A s 



max 



dlnnh 2 



dlnq 



In Figure [10] we evaluated A at two-loop level and presented as a function of the Higgs 
mass, without imposing any restriction on the latter. It turns out that A^ can have acceptable 
values even for close to 120 GeV. In Figure [TT1 the total fine-tuning, defined as max{A, A^} 
is presented as a function of m^. Its value is only slightly larger than that found earlier for 
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Figure 11: Left (Right) figure: The total fine tuning, max{A, A } in function of the Higgs mass, 
at two-loop level, for a la (3a) WMAP bound (in red), respectively. The blue (darker) points sub- 
saturate the dark matter relic density. The continuous line represents the minimal value of the EW 
fine-tuning computed earlier. 

A alone with WMAP saturated dark matter abundance (in red in the plots) . From Figure [11] 
we predict, from minimising max{A, A^} and from consistency with the 3a WMAP bound: 

m h = 114.70 ± 2 GeV, max{A, A n } = 15, (sub-saturating WMAP bound) 

m h = 116.98 ± 2 GeV, max{A, A n } = 28.7, (saturating WMAP bound within 3a) (16) 

The last predicted value of mh is only marginally above that predicted in (fT4"j) . based on 
minimised electroweak fine-tuning and right dark matter abundance. 

5 Predictions for the superpartners from fine tuning limits 

The results so far demonstrate that electroweak fine-tuning has a strong sensitivity to param- 
eters such as fjL, Tni/2, with a preference for lower values. Regarding the mo dependence, A 
has a rather flat dependence when we are in the focus point region. The states that are dom- 
inantly controlled by the fi, m-i/2 parameters are then the most important in determining the 
naturalness of the proposed theory. These include the neutralinos, charginos and the gluino 
states. Further, setting an upper bound on electroweak A gives a bound on the spectrum. If 
any of these states have masses in excess of those given in Table it will require less than 
1% tuning (A > 100) for the MSSM. 
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Table 2: Upper mass limits on superpartners in GeV such that A < 100 remains possible. 

These upper mass limits scale approximately as JAT, so they may be adapted depending 
on how much fine tuning the reader is willing to accept. Overall low fine tuning prefers a 
Higgsino mass of O(0.5 TeV), a gluino of 0(1.5 TeV) and chargino and neutralino masses of 
O(300 GeV). Stop and sbottom masses are significantly larger at 0(3 TeV) due to the weak 
limit on ttiq (focus point). 

Finally we return to the intriguing fact that minimum electroweak fine tuning plus correct 
dark matter abundance corresponds to a Higgs mass just above the LEPII bouncjfl As we 
noted above this point is fixed by the current bounds on the SUSY spectrum and not by the 
current Higgs mass bound which is not included when doing the scans giving Figs [H El [9J [TT1 

One may interpret the SUSY parameters corresponding to this point as being the most 
likely given our present knowledge and so it is of interest to compute the SUSY spectrum for 
this parameter choice as a benchmark for future searches. This is presented in Table [3] where 
it may be seen that it is somewhat non-standard with very heavy squarks and sleptons and 
lighter neutralinos, charginos and gluinos. This has similarities to the SPS2 scenario [36J. 

5.1 Predictions for SUSY searches at the LHC 

It is clear that there is still a wide range of parameters that needs to be explored when testing 

the CMSSM. Will the LHC be able to cover the whole range? To answer this note that, for 

a fine tuning measure A < 100, one must be able to exclude the upper limits of the mass 

parameters appearing in Table [2j Of course the state that affects fine tuning most is the 

8 One may ask whether the fine-tuning measure used above has indeed a clear physical meaning. One can 
object that nature may not choose "minimal" fine-tuning results. One can invoke the example of classical 
chaotic systems, displaying the familiar "butterfly effect" where small variations of the initial conditions bring 
large changes of the final state ("fine tuning"), yet the system is "realised". Such effects exist in (nonlinear) 
classical systems, where initial close values (states) of a parameter exponentially diverge after evolving according 
to the dynamics of the differential eqs. In our setup, one could have such effects not from evolution in time 
but from evolution wrt the energy scale, from the high scale to the low scale, after including quantum effects 
encoded in the RG differential equations. By this analogy one could object that using criteria of low fine-tuning 
to obtain mass bounds (for mh) may not be appropriate. However, the difference is that the discussion in the 
text is at quantum level, so the counterexamples of classical (chaotic) systems do not necessarily apply. 
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Table 3: The favoured Constrained MSSM spectrum of minimal A = 15 giving a sub-saturation of 
the WMAP bound. Masses are given in GeV. 

Higgs scalar and one may see from Figure Q] that establishing the bound > 120 GeV will 
imply that A > 100. However the least fine tuned region corresponds to the lightest Higgs 
consistent with the LEPII bound and this is the region where the LHC searches rely on the 
h — > 77 channel which has a small cross section and will require some 30 fb^ 1 at \fs = 14 TeV 
to explore. Given this it is of interest to consider to what extent the direct SUSY searches will 
probe the low fine tuned regions. Following the discussion in the previous section, the most 
significant processes at the LHC will be those looking for gluinos, winos and neutralinos. 

Studies of SUSY at the LHC [37] have shown that the LHC experiments have a sensitivity 
to gluinos of mass 1.9 TeV for y/s = 10 TeV, 2.4 TeV for y/s = 14 TeV and luminosity 10/6 -1 . 
Of relevance to the first LHC run the limit is 600 GeV for y/s = 10 TeV and luminosity 
100 pb^ 1 . These correspond to probing up to A = 120, 180 and 14 respectively. As we have 
discussed charginos and neutralinos can be quite light, but their signal events are difficult for 
LHC to extract from the background, owing in part to a decreasing — mass gap as 
\fi | decreases [38} I39j . An Atlas study |40j of the trilepton signal from chargino-neutralino 
production found that 30/6 -1 luminosity at 14 TeV is needed for a 3a discovery significance 
for M 2 < 300 GeV and p < 250 GeV [H]. 

6 Summary and Conclusions 

Supersymmetry was introduced to solve the hierarchy problem and to avoid the large fine- 
tuning in the SM Higgs sector associated with the Planck or GUT scale when quantum 
corrections are included. While this hierarchy problem is solved by TeV-scale supersymmetry, 
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the non-observation, so far, of SUSY states means that the MSSM has acquired some residual 
amount of fine-tuning related to unnatural cancellations in the SUSY breaking sector. The 
goal of this paper was to analyse in detail the level of fine tuning in the CMSSM. 

The fine tuning measure A provides a quantitative test of SUSY as a solution to the 
"little" hierarchy problem and measures the "tension" required to satisfy the scalar potential 
minimum condition v 2 ~ —fnlusy I *\ f° r a combination of soft masses m 2 susy ~ TeV, with an 
effective quartic coupling A < 1 and v ~ 0(100) GeV. Although the exact upper limit on the 
fine tuning A beyond which a theory fails to solve the hierarchy problem is debatable, it is 
preferable, for a given model, to have a parameter space configuration corresponding to the 
lowest value of A. We evaluated A at two- loop order and also paid particular attention to 
threshold corrections and to the tan/3 radiative dependence on the parameters. Such effects 
on fine-tuning were not fully considered in the past and turned out to reduce fine tuning 
significantly. 

Our determination of the fine-tuning measure for the CMSSM included the theoretical 
constraints (radiative EWSB, avoiding charge and colour breaking vacua), and also the ex- 
perimental constraints (bounds on superpartner masses, electroweak precision data, 6—7-57, 
b —7- fi [i and muon anomalous magnetic moment, dark matter abundance). As far as we are 
aware, our study is the first two-loop analysis of the fine tuning problem in the CMSSM, 
largely based on SOFTSUSY and micrOMEGAs, SuSpect and our own Mathematica code. 
The latter was very important since it reduced to a feasible level the CPU run time necessary 
to scan the full parameter space. 

Not including the dark matter constraint, we found the minimum value is given by A = 8.8. 
Remarkably, even without imposing the LEP bound on the Higgs mass, the condition fine 
tuning should be a minimum predicts rrih = 114 ± 2 GeV. Adding the constraint on the dark 
matter relic density, one finds A = 15 corresponding to rrih = 114.7 ± 2 GeV and this rises to 
A = 17.8 (m h = 115.9 ± 2 GeV) for SUSY dark matter abundance within 3a of the WMAP 
constraint. The results are encouraging for the search for SUSY because we considered the 
"conservative" case of CMSSM, and it is well-known that relaxing gaugino universality can 
reduce A further [6j |9] . 

The spectrum corresponding to the minimum value of the fine tuning shows similarities 
to the SPS2 scenario with light neutralinos, charginos and gluinos (corresponding to light 
£*> m i/2) an d heavy squarks and sleptons corresponding to large mo, near the focus point 
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limiting value [13} I14j . It provides the "best" estimate for the SUSY spectrum given the 
present experimental bounds. 

Increasing nih above the minimum fine tuned value causes A to increases exponentially 
fast and one leaves the focus point region at the edge of which this minimal value is reached; 
one obtains A = 100 (1000) for a scalar mass rrih = 121 (126) GeV, respectively. Ultimately 
the question whether the SUSY solution to the hierarchy problem has been experimentally 
tested relies on what value of fine tuning represents the limit of acceptability. Given a value 
one can determine the range of parameter space that is still acceptable. For the case that the 
fine tuning measure should satisfy A < 100, we determined the corresponding superpartners 
masses and CMSSM parameters values, that can be relevant for SUSY searches. 
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7 Appendix 

A.l Higgs mass and EW fine tuning 

We provide technical results used in the text to evaluate A. The potential used in (P) 



V 



m\ \Hi\ 2 + ml \H 2 \ 2 - (m| #1 • H 2 + h.c.) 



+ n -Vl |#l| 4 + - A2 1 4 + A3 |#l| 2 \H2\ 2 + A4 \H\ ■ H 2 \ 2 



+ 



1 



A 5 (#1 ■ H 2 ) 2 + A 6 |J^i| 2 (#1 ■ H 2 ) + A 7 \H 2 \ 2 (# x ■ H 2 ) + h.c. 



(A-l) 



with H1..H2 = #1 #2 - H i H 2 - Using the notation 



m 



rri! cos /3 + m 2 sin /3 — m 3 sin 2/3 
Ai 



A, 



A = y cos 4 /? + y sin 4 /3 + ^ sin 2 2/3 + sin 2/3 (A 6 cos 2 /3 + A 7 sin 2 0) (A-2) 



the minimisation conditions give 

v 2 = —m 2 /X, 

or, equivalently, 



9m 2 2 <9A 

2X W = m w 



(A-3) 



; - ™i + ™2 + ^- [A1C 2 +A 2 s| + A 3 45 + (A 6 + A 7 )s 2( g + A 6 cot/3 + A 7 tan/3] 



sin 2/3 



m 2 — m 2 tan 2 /3 = — — [cos 2 /3 (A x — A 2 tan 4 /3) + sin 2/3(A 6 -A 7 tan 2 /3)] (A-4) 



One finds for the CP odd Higgs mass: 



2/7? ip 1 

— -h-J (2A 5 + A 6 cot/3 + A 7 tan/3) 
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(A-5) 



with g 2 = g 2 + g 2 - In the notation of [42] . the CP even Higgs masses are 
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(A-6) 



where A 345 j = A 345 — A,-, A is that of ([2]), and, (sp,cp) = (sin /3, cos /3). In the limit of large 
m A , m 2 reaches an upper limit of 2Xv 2 (which tends to A 2 t» 2 for large tan/3). At tree level, 
A = (g 2 /8) cos 2 2/3, and the individual A„- are: 

,2 



^1,2 — — -^345 — ^5,6,7 







(A-7) 



The general formula for fine-tuning A p wrt a parameter p in a two-higgs doublet model that 
we are using can be found in the Appendix of [IT] . For large tan /3 = V2 jv \ this reduces to 
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<91ni; 2 
<91np 
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3452 I 



where x' = dx/dlnp is the partial derivative of x wrt p. 



A. 2 The scalar potential: 1 Loop Leading Log (1LLL) Terms 

Here we show how to obtain the one-loop corrected potential, which is "improved" to the 
two-loop leading log (2LLL) result in the next section. Start with 

V® = m 2 ^ 2 + fh 2 \H 2 \ 2 -fh 2 {H 1 H 2 + Yi. C .) + 9 ^- ({H^ 2 - \H 2 \ 2 Y (A-9) 



This receives (field dependent threshold) corrections, computed using the Coleman- Weinberg 
potential [33]: 

y {1) = ^2E(- 1 ) 2Jfe ( 2 ^+ 1 )^^(iog^|-0 (A-10 

where m k is the field dependent mass, the degeneracy factor g k is 6 for squarks, and J k is the 
particle spin. All parameters in eq, (|A-10p are evaluated at the scale Q using the RGEs which 
ignore the particle masses. The field dependent squark masses are (neglecting 0(g 4 ) terms): 
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(A-ll) 
(A-12) 



and where mg ur>{M s ) = M s is assumed. Here, M s is the soft SUSY breaking squark mass 
evaluated at the squark mass scale. 

One can expand the non-linear field dependence (log) in in inverse powers of I/-M5 
to find the dominant threshold corrections, which come from the third generation squarks: 
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t = log(M 2 s /Q 2 ' 



(A-15) 



When running below the EWSB scale, the inclusion of higher dimensional terms (threshold 
corrections) lead to a re-summation such that M s is replaced by a mass scale related to the 
physical particle masses |44| . For the results of this paper, the geometric mean of the stop 
masses is used in the place of Ma. 

The above equations are valid down to the top mass scale; below this scale threshold 
corrections from the top quark should also be included. The dominant effect of running below 
the top scale can be absorbed by setting Q in the above equations as the "running" top mass 
evaluated at the scale Q instead of the pole mass. 
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From eqs. (|A-9j) to (|A-14p one obtains the parameters entering in the scalar potential ([I]), 
evaluated at the scale Q (below Ms), in the one-loop leading log approximation (1LLL): 
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3/> 4 / /x 2 _ /A4|\ 3hi (» 2 _ /A4 2 \ 

T ic_2 »/r2 o«4 J "I" ic_2 71^2 otit4 ^ Zi J 



16vr 2 \M 2 3M|7 16vr 2 \M 2 3M| 



3/> 4 n 2 A 2 3ht v 2 A 2 



96?r 2 M| 96vr 2 M 4 



(A-22) 



g 2 / 3^ (h 2 A b - h 2 A t ) \ 3hj ^A t 3hj p f A* 6A b 

6 4^ 32vr 2 M 2 J 96ir 2 Mj 967r 2 M s \Mg M s } ' 

g 2 ( 3»{h 2 t A t -h 2 A b ) \ 3hj y?A h 3hj » (A\ 6A t 

7 4^ 32vr 2 M 2 J 96tt 2 jVfJ 967r 2 M s m s. 



These analytic results agree with [25J which ignore the stop mixing corrections to the D-terms, 
but are included here for completeness. The following notation is used in this appendix. 

2A tb ( A tt \ 
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A. 3 The scalar potential: 2 Loop Leading Log (2LLL) Terms 

The two-loop leading log (2LLL) Coleman- Weinberg potential can be found in the arXiv 
version of 05] to O (g% hf , g% hf) and 0(h$,hf h\,h\ h%,h 6 b ) , see also [Ml 07] for the general 
case. The method of the previous section may be used to determine the 2LLL contributions 
to the Higgs scalar potential, however here we proceed instead with an approach similar to 
that in [44J, to RG-improve the 1-loop result into a 2LLL result. A step approximation is 
applied to the /3-functions so that the MSSM RG eqs are used between the GUT and stop 
mass scale, then the 2HDM SM RG eqs between the stop and top mass scales, and finally the 
top is integrated out to reach the electroweak scale. 

When setting the renormalisation scale in eqs (|A-16|) to (jA-24|) as Q = M s , the logarithmic 
terms are removed but the finite corrections from stop mixing remain. These results are then 
used as boundary conditions for the parameters at the scale M s . A series expansion of the 
RG eqs is then applied: 



the Coleman- Weinberg potential approach. For a /3 A -function of the form bX + c, eq (|A-27j) 
becomes 



where the couplings are evaluated at the scale Q unless stated otherwise. The /3-functions for 





(A-28) 



the 2HDM SM |4]8] are listed below, neglecting 0(/i^) terms, and with the (3 X -functions also 
neglecting 0(<? 4 , g 2 X i , A?) terms: 






I {hi + hi) ml + O ( 9 2 m 2 ) 



(A-29) 
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a 6/) 2 fA, - hi) 


16vr 2 /3 A2 p 


a 6/ i 2 (A 2 -/ tt 2 ) 


16tt 2 /3 A3 p 


a 3A 3 (/ l 2 + / i 2 )-6/ l 2 / l 2 


16vr 2 /3 A4 p 


a 3\(h 2 + h 2 b )+Qh 2 h 2 


16vr 2 /3 A5 * 


a 3A 5 (/ l f + ^) 


16vr 2 /3 A6 p 




16vr 2 /3 A7 p 





and finally 



16vr 2 ^ 



IGtt 2 /?,, 



*2 



2^+2^ 



8^1 



f2 



17 



/* 2 + ^ 2 + ^ 2 



8^1 



91 



9 2 5 



12 

■al 



9i 



(A-30) 

(A-31) 
(A-32) 



2 2 1 T ° 4 z 12 

Using ()A-28p . the analytic 2- loop results in [25] are then recovered when the same level of 
approximation is considered. For example, 

3h 2 



An 



[A 2 (M s ) -X 2 a 2 t]-b 2 t + 

3hf 



X 2 (M s )-X 2 ^j-t 



+ 



8vr 2 



a 2 6 2 + 
t + 



16vr 2 
1 



2^ ftf " /3a 2 ) + 0(A) 



16vr 2 



5« + 5« 



(A-33) 



(A-34) 



The couplings entering in the expression of A 2 (Mo) are re-expressed in terms of their values 
at the scale Q, (with a logarithmic correction which compensates for the running below M s ): 



h$(M s ) 



hUi 



16vr 2 



(9h 2 + h 2 -16g 2 )+0{g 2 t,t 2 ) 



(A-35) 



hf(M s ) = hi^ + ^(Vh 2 + h 2 -lQg 2 )+0{g 2 t,t 2 )} (A-36) 



This leads to the following expression, in agreement with [25J, when the stop mixing contri- 
butions to the D-terms in the potential are neglected: 



A 9 



9 



1 %t 

8vr 2 



96vr 2 M| 



t 



1 + ^ {9h 2 b - 5h 2 - 16g. 



+ 



3hf 
8tt 2 



t+ 2 + 16vr 2 



t (3h 2 hi 



+ ^-8g 2 )(X t + t) 



(A-37) 
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Note that these results assume that the CP odd Higgs mass is not decoupled. If this is the 
case, the usual SM /3-functions should be used. The effective quartic coupling at the EW scale 
when m A < M s is given by: 



A « ^ cos 2 2/3 [l-^ (h 2 b +h 2 t + (h 2 b -h 2 t )sec2/3)t 



3 h 4 

° 'H -An 

1 h 4 

+ 16^ C ° S P 



X, 



1 



t + ~2~ + IGtt 2 



16vr 2 
1 /3/i 2 



3A? 
2 2 



^ 2 

h — 

2 2 



8^3 



8 5 |) (x fc t + t 2 ) +( 5 2 



(A-38) 



with the following notation: 



<5i = 



3t(h 2 b -h 2 ) i t /xcot/3 



16vr 2 



3t(h 2 -h 2 ) A h n tan /3 



16vr 2 



6M 2 



(A-39) 
(A-40) 



where X tb is defined as X tb (A fb -> I t 6 ) with 



A t — n cot /3 
A b — n tan /3 



(A-41) 



A similar but distinct result is obtained when ~ M s (notably no S i terms and a differ- 
ent dependence on tan/3 and the mixed Yukawa couplings). The threshold corrections are 
dependent on where the CP odd Higgs decouples. The same procedure has been applied to 
determine the 2LLL threshold corrections to the mass terms. 
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